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Abstract 

We derive expansions of the Hermite and Laguerre kernels at the edge of the spectrum of 
the finite n Gaussian Unitary Ensemble (GUE„) and the finite n Laguerre Unitary Ensem- 
ble (LUE„), respectively. Using these large n kernel expansions, we prove an Edgeworth 
type theorem for the largest eigenvalue distribution function of GUE„ and LUE„. In our 
Edgeworth expansion, the correction terms are expressed in terms of the same Painleve II 
function appearing in the leading term, i.e. in the Tracy- Widom distribution. We conclude 
with a brief discussion of the universality of these results. 

1 Introduction 

The limiting distribution function for the largest eigenvalues in orthogonal, unitary 
and symplectic random matrix ensembles have found many applications outside their 
initial discovery in random matrix theory, see, for example, [HI IH 123 recent re- 
views. In these applications it is important to have correction terms to the limiting 
distribution. For example, in statistics jT2] the sample size is always finite; and to 
assess quantitatively the range of validity of limit laws, one needs finite n correction 
terms. In classical probability, a similar issue arises in the application of the Cen- 
tral Limit Theorem (CLT) to finite n problems. Here the two main results are the 
Berry-Esseen theorem and the Edgeworth expansion jH]. 

Recall if Sn is a sum of i.i.d. random variables Xj, each with mean /i and variance 
0"^, that the distribution F„ of the normalized random variable {Sn — nfi)/{ay/n) 
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satisfies the Edgeworth expansion^ 



= (/>(x)^n-5^+^i?j(a:) + o(n-^"+^) (1.1) 
i=3 

uniformly in x. Here $ is tfie standard normal distribution with density 0, and Rj 
are polynomials depending only on E(Xj^) but not on n and r (or the underlying 
distribution of the Xj). 
Introduce 

Fnh\t)=^GA>^^L<t) (1-2) 

where A^;^ is the largest eigenvalue in GUE„ or LUE„, respectively. (When the 
context is clear, we often drop the G or L.) To obtain a nontrivial limit theorem, we 
must, as is well known, define normalized random variables We find it useful 

to "fine tune" our normalization (see also 



Xg,^ - {2{n + CG)f' 
,L ^ ^max - 4(n + cl) - 2a 

^max • 2(2n)V3 ^ • ) 

where cg,l are constants (to be chosen later), and a is the parameter appearing in 
LUE„. (That is, the parameter a appearing in the Laguerre polynomials L".) Then 
•^m!M converge in distribution to GUE Tracy- Widom (commonly denoted F2). In this 
paper we initiate the study of Edgeworth expansions for both GUE„ and LUE„; that 
is, we find the analogue of (jl.lj) for F^'^ — F2. We now state our main results. 

Our first result is an extension of the Plancherel-Rotach theorem for the Laguerre 
poynomials We set 

1 X 

^ = {in + 2q; + 2c)2 H — 2 — r where X and c are bounded, (1.5) 

23?T,6 

and denote by Ai the Airy function (see, e.g., jlHI)- 
Theorem 1.1. For a > — 1 we have as n ^ 00 

e-^'/^L'^ii^) = (-l)"2-°-^n-^ ( Ai(X) + ^^-^ Ai'(X)n-^ + 

I 23 



X Ai(X) H 2- Ai'(X) 



10 -21 ' ' 20-2i 

,5a — 15ca + 2c^ — 15c^ — 56c — 6 , c — 1 



_ 2 
n 3 + 



(c - 1) (5(c - 2)c - 3(2 + 5a)) ^ 



60 



60 40 

+ 0(n-t)Ai(X) j 



X') ki{X) + 



^Wc assume, of course, the moments E(Xj^), k = 3,...,r, exist; and as well, the condition 
lim|(^l^oo sup |v(C)l < 00 where Lp is the characteristic function of Xj, see [H]. 
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From Theorem 11.11 we derive expansions for both the Hermite and Laguerre ker- 
nels. Recall that if^ 



HJx) e-"'/2 and 6'^(x) = x"/' e""/' L°( 



then the Hermite kernel is 

n-l 



Kn{x, y) = Vk{x)(pk{y) (1-6) 



k=0 

and the Laguerre kernel is 

n-l 



A;=0 

Finally, the Airy kernel is 



K'^(xv) = T '^^^^^'^^^^^ (17) 



, , Ai(x) Ai (y) - AUy) Ai (x) ^ , , , , , , , , 

KAi{x,y) = ^^ \^ =1 Ai{x + z)A:i{y + z)dz. (1.8) 

X ~ y Jo 

Using Theorem II. II we prove 

Theorem 1.2. For x = (2(n + cg))^ + 2"5n~ix and y = {2{n + cg))^ + 2-^n~W 
with X, Y and cq bounded, 

K^{x, y) dx = [Kp,^X, Y) - CG Ai{X) Ai{Y)n-h + 
^ [{X + Y) Ai'(X) Ai'(F) - (X^ + XF + Y^) Ai(X) Ai(F) + 



-Orir^ -I- S 

— ^— (Ai'(X) Ai(F) + Ai(X) Ai'(r)) 



n~i + 0{n~^)E{X, Y) | dX. (1.9) 



The error term, E{X, Y), is the kernel of an integral operator on L'^{J) which is trace 
class for any Borel subset J of the reals that is bounded away from minus infinity. 

For the Laguerre kernel we prove 

Theorem 1.3. For x = A{n + cl) + 2a + 2{2n)^ X and y = 4:{n + cl) +2a + 2{2n)^iY 
with X , Y and cl bounded, 

K{x,y) dx= |j^Ai(X,r) -2icLAi(X)Ai(F)n-3 + 



^Here Hn{x) are the Hermite polynomials of degree n. 
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^ [ (X^ + XY + Y^) Ai{X) Ai(Y) -{X + Y) Ai'(X) Ai'(F)- 

(lOci - l)(Ai(X) Ai'(y) + Ai'(X) Ai(r))] n-i + 0{n-^)F{X, F)} dX. (1.10) 

The error term F{X, Y) is the kernel of an integral operator on L'^{J) which is trace 
class for any Borel subset J of the reals which is bounded away from minus infinity. 

To state our main theorem, we need a number of definitions. First define the 
constants 



acG,2 = cg, <,2 = 2^Ci, = b^ = ^. 

Following the notations of Tracy and Widom [201, 



^cgA-^) = - 3u2 + {-20cq + 3)vo + uiVq - uqVi + uqVq - uIwq, 
^clA-^) = - 3u2 + (20c^ - 2)vo + uiVq - uqVi + uqvI - uIwq, 

where 

/»oo poo 

Ui := Ui{s) = I q{x)x^ Ai{x) dx, Vi := Vi{s) = / q{x)x'^ Ai{x) dx and 

J s J s 

POO 

Wi := Wi{s) = / q'{x)x'' Ai'{x) dx + Uo{s)vi{s) 



with g the solution of the Painleve II equation, q" = sq + 2q^, subject to the boundary 
condition q{s) ~ Ai(s) as s — > oo. Finally, the GUE Tracy- Widom distribution is 

F2(s) = det(/-i^AiX(,,oo)) =exp (^-^ {x - s)q{x f dx^ (1.11) 

where K^i is the operator with Airy kernel (jl.8p and X(s,oo) is the indicator function 
of the interval (s, oo). With now can state our main result. 

Theorem 1.4. We set 

t= {2{n + CG))^ +2-^n-^ s for GUE^ (1.12) 

and 

t = 4:{n + cl) + 2a + 2{2n)^ s for LUE^. (1.13) 

Then as n oo 

F^tit) = F2{s){l + a5^^,2 ^o(^) + b^'^^E^'^^.is) n-l} + Oin~') (1.14) 
uniformly in s. If in addition, 

4 + 4 = ^) then E^^ ^is) = £'f;,2(s) = ^c,2(s), and 
Fn,2\t) = F2is){l + a5^^,2 uo{s) n--^ + 6^'^i5c,2(s) + 0{n~^). (1.15) 
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Note that the n~2 correction term in the Edgeworth expansion for the CLT is 
universal in the sense that only the constant factor depends on the underlying distri- 
bution. We see in (jl.isp a similar universality, and we conjecture that this universality 
extends to a wider class of unitary ensembles. 

In §2 we derive Theorems 11.21 and 11.31 from Theorem 11.11 In §3 we follow [20] 
to prove Theorem 11.41 The proof of Theorem 11.11 will be given in the Appendices 
together with some facts needed to prove the last Theorem. 

2 Correction terms for the Hermite and Laguerre 
kernel at the edge of the spectrum 

To simplify notations we will use matrix ensembles of (n + 1) x (n + 1) matrices 
throughout this section and part of the next section. After the fine tuning of the 
variables in ^3.31 we will use ensembles of n x n matrices. 



2.1 Hermite case 

We have the following representation of the Hermite kernel from the Christoffel- 
Darboux formula. 



Kn+i{x,y) = ^^- ^-^ (2.1) 

As mentioned in SJQ Theorem II. 21 is a corollary of Theorem ll.il We recall the relation 
between the Hermite and Laguerre polynomials 

H2n{x) = (-l)"22"n!L;^(x2), H2n+i{x) = (-l)"22"+in!xL|(x2) 

so as to estimate the right side of 1)2.11) . We will also assume^ without lost of general- 
ities that n = 2k is even. Using the symmetry of the Hermite kernel in 1)2.1)) we only 
need to find an expansion of 

,. 2''^+\k\fxe-4Ll{x^)e-4Lf'{y^) 

u:>n+i(x)(pJy) = — — — 1 (2.2) 

v.„+n J^nKyj (24'=+i(2A; + l)!(2fc)!7r)^ ^ ^ 

since the other term follows by interchanging x and y. The Laguerre polynomial of 
argument y in ()2.2)) has parameter a = — |; thus, in order to apply Theorem 11.1) we 

write y in the form y = (Ak — 1 + 2c')^ -\ — ^-t which corresponds to c' = c -|- 1. Next 

we use Stirling's formula to estimate 

^ ^ -(-l)^2"2-iA;"3(_i)'^25-3A;-3 . 



(24fc+i(2A; + l)!(2A;)!7r) 
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For n odd the same analysis can be carried out to produce the same result. 
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The last factors in the left hand side are the constant factor in Theorem (ll.lj) for the 
Laguerre functions in x and y respectively. This expression is 



^ 2ikl 4:k 



This times, the constant^ from (j2.ip . gives 



2U^(l + ^ + ^ + 0(A;-t)). 



23 A;3 



Ak 



We substitute all these into ()2.2p . and then interchange x and y to have the second 
term ipri+i{y)^n{x). Finally with the help of Mathematica we derive the following 
version of Theorem 11.21 



For X = (2n + 1 + 2c)^ + -^-^ and y = {2n + 1 + 2c)2 + -^-^ , (2.3) 

22n6 22?T,6 



i^„+i(x,i/)dx= <^iCAi(X,F) 



1 - 2c 



Ai(X) Ai(F)n" 



X + F 
20 



Ai'(X) Ai'(r)- 



+ + AKX) AKF) + ^1^^^±^(A.'(X) MiY) + AKX) A.'(y))' 



20 V / V / 20 



_ 2 

n 3 



,1 



+ 0(-) E(X,F) UX (2.4) 
n ) 

Using the symmetry between the X and Y terms from ^2.1^ . this error term can be 
expressed finite sum 

E{X,Y) = ^P,(X,r)Ai(X)Ai(r) + 5^Q,,(X,F)Ai'(X)Ai(F) + 

j jl 

J2 Qn y) Ai(^) Ai' iY) + Ai' {X) Ai' (F) + 

32 j 

X^Y^ Ai(X) Ai' (F) - X^F-'' Ai' (X) Ai(r ) 



x-r 

where all coefficients of the polynomials P, Q and i?, and the aj^k have a factor of 
A: e {0,1,2,3}. 

The first four terms are kernels of finite rank operators on any Borel subset J of 
M not including minus infinity;^ and thus, are trace class. The last term also defines 

'''Recall that in this section the notation Kn stands for an ensemble of (n+ 1) x (n + 1) matrices. 
^This last restriction is due to the behavior of the Airy function near minus infinity. 
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a trace class operator. This is best seen from the following result. If we assumed 
without lost of generalities that j < k, and set k — j = s, then 

X^Y'' Ai(X) Ai' (Y) - X^Y^ Ai' {X) ki{Y) 
X-Y 

{Y^ -X^ + X-) Ai(X) Ai'(y) - X- Ai'(X) Ai(y) 

= x^r 

/ s-l ^ 

= {XYy {x'Km{X,Y) + ^y''-*X^Ai(X) Ai'(F) 
This shows that the error E{X, Y) in ()2.4|1 is the kernel of a trace class operator. 



2.2 Laguerre case 

Again by the Christoffel-Darboux formula, 

^-+1^^' y> - ^ n + a\ ' ^ ' 



Via 



n 



2.2.1 Asymptotic of e"iL^(x) at x = 4n + 2a + 2c + 2(2n)^X 

In order to apply Theorem 11.11 to the Laguerre kernel at the edge of the spectrum 
(corresponding to a; = 4?7, + 2a + 2c + 2(2?7,)3X for bounded c and X), we need to 
make a variable change x = where ^ = A/4n + 2a + 2c + 2~3?t,~6 t. We use 
these two expressions to solve for t in terms of X and then substitute this value for t 
into Theorem 11.11 to obtain the desired asymptotics. But^ in order to have accurate 
asymptotics for the Laguerre functions at the edge of the spectrum, we will use the 
expression of ^ involving = (4n + 2a + 2c) Thus 

x = i^ is equivalent to + 2(2n)^X = + (2/„)it + (2/„)"§t2 . 



This quadratic has solutions t± = -(-(2/„)t ± (2/„)i y 1 + 4(2/„)"22(2ra)iX) . 

Actually, only the solution with the plus sign is to be taken as it is the only one 
bounded when n increases. An expansion of = t leads to: 

(a + c) X^ ^ , 4, 
t = X- ^-r^X - + 0(n-3) . 

6n 23723 

Thus if in Theorem II. II we replace X by this value of t, we obtain, again with the 
help of Mathematica, the desired expansion. 

^See the last footnote in the proof of Theorem ll.ll m the Appendix for this technical point. 
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Lemma 2.1. For x = An + 2a + 2c + 2{2n)3X and X bounded, 



e-iL^(x) = (-l)"2-"-3n-3 } Ai(X) 



23 



—XAi{X) jAi (X) 



10 ■ 23 



10 ■ 23 



n 34- 



— [- (6 + 56c + ISc^ - 2c^ + 5a(3c - 1) - GX^ + QcX^) AUX) 
60 



(6 + a(5- 15c) -6c- 15c2 + 5c^)XAi'(X) + 0(n"^/^) Ai(X)| (2.6) 



2.2.2 Asymptotic of e"tL^+i(x) at x = 4n + 2a + 2c + 2(2n)^X 
Making use of this last formula, we can derive an asymptotic for 

X 1 

e~2L°_,_;^(x) when x = 4n + 2a + 2c + 2(2n)3X. Note that the degree of the Laguerre 
polynomial is no longer n but n + 1 , so in order to use Lemmg lTiTl we need to write x 
in terms of n + 1 or x = 4(n + 1) + 2a + 2(c - 2) + 2(2(n + 1))^ [X - ^X + 0{n-^)]. 
The substitution needed here is c ^ c — 2 , and X — > X — -^X. Matheniatica again 
gives: 

Lemma 2.2. For x = An + 2a + 2c + 2(2n)^X and X bounded, 



-l)"+^2-"-t(n + l)'t \ Ai(X) + i^-^ Ai'(X)n-i + 

2 



42 - 30c + 5c2-5a^ ^ 2X2 ■ 
— XAi(X) 2-Ai (X) 



10 ■ 23 



10 ■ 23 



n 3 + 



60 



[(30 + 28c - 27c2 + 2c^ - 5a(3c - 7) + 18X=^ - 6cX3) Ai(X) 



+ (-102-5a(3c-7) + 114c-45c2 + 5c3)XAi'(X) + 0(^-^/3) Ai(X)} 



To complete this subsection we need to estimate 

{n + l){xy)^ 



rra + 11 



n + a 
n 



We have 



(xy)2 = 22"n"[l + 



a{X + Y) a(a + c 



+ 



2tni ' 2n 



(2.7) 



0(n-3)], 



n + a \ a'^ + a 



r(a + l)^ - j=n'^[l + ^— + Oin-^)] 
and in addition, the product of the constant factor on the right of ()2.6|1 and ()2.7p is 



;-l)"+^2-"-3(n + 1)^3 . (-l)«2-"-3n-3 = -2"'"n^3[i + 0{n-')] 



thus 



r(a + 1 



n + a 



n 



2 1, aiX + Y) 3a(c-l)+4 4,, 

-2-3^3 1 + ^ ^ , > + ^ ^ + O n-3 

23/13 6r2 

Substituting all these quantities in ()2.5|) . give the following version of Theorem 11.31 
For X = 4n + 2a + 2c + 2(2n)^X and y= 4n + 2a + 2c + 2(2n)5r (2.8; 

K;^+i(a;, y)dx = {km{X, Y) + ^ Ai(X) Ai(F) + 

I 23 



2{X^ + XY + Y^) Ai(X) Ai(r) - 2{X + Y) Ai'(X) Ai'(F)- 



2i 10 



(18-20c + 5c2)(Aim Ai'm + Ai'm Ai(F))l n-i + Oi-)F{X,Y)\dX. 

\ n J 

(2.9) 

As in the Hermite case, the error term F{X, Y) is a finite sum of kernels of trace class 
operators, therefore is a kernel of a trace class operator on L'^{J) for any subset J of 
the reals which is bounded away from minus infinity. 

2.3 Conclusion 

For our order of expansion, we see that both kernels are finite rank perturbation of 
the Airy kernel. In the Laguerre case, the final result does not involve an explicit 
presence of the order a. 

3 Expansion of the Fredholm determinants at the 
edge of the spectrum 

This part will only make use ()2.3|) and ()2.4j) to derive the desired result in the Hermite 
case and ()2.8|) and ()2.9|) in the Laguerre case. Most of the derivations will follow from 
the work of Tracy and Widom. 
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3.1 Edgeworth expansion of 

Recall that 

F^^.^^it) = det{I - Kn+i) (3.1) 

where this determinant is the Fredholm determinant of the integral operator with 
kernel ii'„+i(x, y) on L^(s, oo), t and s are related by 

t = (2ri+l + 2c)^ + 2~^n-^s. (3.2) 

In this section we will estimate this determinant. Most of our derivations involve 
trace class operators where the Fredholm determinant is analytic. If in p.ip we use 
the expression of -R'n+i given by (j2.4|) . the continuity of the determinant (in trace 
class norm) allows us to pull the error term involving the kernel -E(X, Y) out of the 
determinant as an 0{n~^) term. We therefore have 

det(J - Kn+i) = det ( / - \ KAiiX,Y) + Ai(X) Ai(F)n"U 



^ + ^ Ai'(X) Ai'(F) - + + Ai(X) Ai(F)- 



_ 20 V ^ V / 20 
-lOc^ + 10c- 1 



-(Ai (X)Ai(F)+Ai(X)Ai {¥)) 



n~H)+0{n~^). (3.3) 



20 

Note we are using the obvious notation of writing the kernel for the operator appearing 
in the determinant. We continue to employ this notation below. If we factor out 
det (J - KAi{X, Y)) and set 

Kn+i{x, y) dx - KAi{X, Y)dX := L{X, Y) dX + 0{n~^)E{X, Y) dX, 

the multiplicative property of the determinant gives 

det(/-K„+i(x, y)) = det(/-i^Ai(X, Y))det{I -{I - K^iiX, Y))~^L{X, Y)) + 0{n-^). 

(3.4) 

The first factor on the right is known, so we need to estimate the second factor. We 
will express this factor as a finite sum of rank one operators 

k 

(J - J^Ai(X, Y))-^L{X, Y) = Y, Ol^{XmY) (3.5) 

4=1 

and use the well known formula 

det f/-^a,(X)A(F)j = det (5,, - (a„ /5,))^<^_ (3.6) 
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to transform the problem to one involving estimations of inner products. Let (■, ■) 
denote the inner product in L^(s, oo). A good reference for much of the results that 
follows is [20] • We have the following representation. 

(J - KAiiX, Y)r^L{X, Y) = ^^(/ - Km{X, Ai(X) Ai(F)n-i + 



^ [(/ - Km{X, Y))-^X Ai'(X) Ai'(r) + (/ - Km{X, Y))-^Y Ai'(X) Ai'(r)- 

(/ - Km{X, Y)r^X^ Ai(X) Ai(F) - (/ - Km{X, Y))-^XY ki{X) Ai(F)- 
(/ - AAi(X, Y)Y^Y^ Ai(X) Ai(F) + (-10c' + 10c - 1)(/ - AAi(X, F))"^ Ai'(X) Ai(F) 



+ (-10c' + 10c - 1)(/ - A:Ai(A:, Y))-^ Ai(X) Ai'(F) ^"3 + ©(n"^) 

In the notations above, we think of all the quantities involved as kernels of integral 
operators, the analogues of those in [201 • We therefore have: 



^^(/ - A'Ai(X, Ai(X) Ai(F)n-^ = ^^Q(X) M{Y)n-l := a^{X)P,{Y) 

where Q{X) is the action of the integral operator with kernel (/ — A'Ai)^^ on Ai(X). 
In the same way we have, 

^(J-i^Ai(X,F))-^XAi'(X)Ai'(F)n-i = ^{{I - K^{)-' f){X) k\ {Y)n-i 

:= a2{X)(32{Y) 

where f{Z) = Z Ai{Z) and, 

^{I-KA;iX,Y)r'Ai{X)YAiiY)n-l = ^P(X)/(F)n-i 

:= as{X)Ps{Y) 

where P{X) is the action of the integral operator with kernel (/ — A'ai)"^ acting on 
Ai'(X). 



-^{I-KUX,Y))-'X'Ai{X)Ai{Y)n-l = -!((/- A'^„.,)-i^?)(x) Ai(y)n-i 

:= a,{X)(],{Y) 



where g{Z) = Z^ Ai{Z) , 

-^(/- AAi(X,F))"iXrAi(X)Ai(F)n-t = -^((I - KAr'h)iX)hiY)n-l 

:= a5{X)f],{Y) 
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where h{Z) = Z Ai{Z) , 

-^(J - Ka,{X, Y)r' Ai{X)Y' Ai(F)n-i = -l-Q^X)9{Y)n~l := ae{X)(3,{Y) 

C{I- KAi{X, Y))-^ Ai'(X) Ai(r)n-§ = CP{X) Ai(r)n-5 := a7{X)(37{Y) , and 

C {I - KAi{X,Y))-' Ai{X) Ai {Y)n-I = CQ{X)Ai{Y)n-"^ := as{X)Ps{Y) 

where C = (-lOc^ + 10c - l)/20. 

If we set {ai,f3j) = aijn~s for j ^ I and = aun' 3 , expanding with 

respect to the first row leads to the folfowing expression. 



(1- 




1 -2 
1 - a22n 3 

2 

-032^ 3 


2 

-023"'^ 3 
1 - 033^-3 ■ ■ 


2 

-a28n 3 

-038^-3 






-as2n~^ 


-a83n-i ■ ■ 


■ 1 - a88n-§ 


8 

E( 


-l)Vfc?^ 


-idet(Ci,--- 


5 Ck, ■ ■ ■ , Cg) 




k=2 










(1- 




1 - 022"-"^ 
2 

-032?! 3 


-a23n~^ 
1 - a33'T'"3 ■ ■ 


-a28n"i 
-a38n"3 






2 

-082^ 3 


-a83n"i 


■ 1 - assn'^ 



We factor out from column Ci in the last step. The determinant in the last line 
is of the same form as the the original determinant, therefore a similar transformation 
to this last determinant leads to the following result. 



k=2 



det {6ij - {ai,(3j))^^. .^g = {l-aiin 3)]^(l-afcfcn ^ + 0{n ^) 

*)E(-i 



l-a„n-t)2^(-l)'„-- J2 Yliify+Oin-^) 

k=0 ,ifc;j^^ise{2,.. .,8} j=l 

8 8 

1 - aiin"^ - n-i akk + 0(n"^) = 1 - ^(a^, pk) + 0{n-^). 



k=2 



k=l 



Thus we only need to compute the inner products (at, Pk) for /c = 1, ■ ■ ■ , 8. 

To simplify notations we will write for example uq instead of uo{s). We therefore 

have: 

(q!i,/?i) = — - — {Q,Ai)n 3 = — - — uqU 3 
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(«2,/?2) = ^((/-i^Ai)-^XAi',Ai') = ^(P,XAi') = ^^, 



_ 2 _ 2 

T? 3 / T) 3 



71 3 T? 3 77 3 

(«4,/34) = -^((/-AA0"^^7,Ai) = -^(g,X2Ai) = -^n, 



77 3 77 3 

(«5,/35) = -^((/-AAi)-^/^,/^) = -^(Qi,XAi) 



(a6,/36) = -^(g,X2Ai) = -^M2 



(a^,/?^) = Cn-3(P,Ai) = Cn-3(Q,Ai') = Cn~-^vo 
To estimate (asj/Ss), we use equation 2.12 of jSOj which says 

Qi{X) = XQ{X) + uoP{X) - voQ{X), to have 

iQ^{X),XAi{X)) = {XQ{X)+uoP{X)-voQ{X),XAi{X)) 
= U2-voUi + uoiPiX),XAi{X)) 

= U2- VqUi + UqVi 

= U2 — VqUi + UqVi — UqVq + uIwq. ( We used (2.14) of (201) 

And 

_ 2 

("5, Ps) = ^('"2 - ^l^'O + UqVi - UqVq + ulwo). 

Substituting this into the formula for the determinant gives 

8 

det{Sij - (ai, Pj))i<ij<8 = 1 - ^{ak,Pk) + 0{n~^) 

k=l 

_ 2 

1 — 2c 1 ra 3 , „ 
= 1 —uon'^i - {2wi - 3u2 + (-20c^ + 20c - 2)vo + 

- uqVi + UoVq - ulwo] + 0{n~^). 

If we set 



C / \ ^ 1 C 1 

a^n s = , 09= and 

c,2V ; 2 ' 2 20 



E^2{^) = 2wi — 3u2 + (— 20c^ + 20c — 2)vo + uiVo — uqVi + uqVq — uIwq, (3.7) 
equation (j3.4p gives 



13 



Lemma 3.1. For = {2n + 1 + 2c) 2+2 ^ n e $, 

Fn+iA^c) = F^is) {l+a%Uo{s)n-^+b^E^^,is)n~l]+0{n-') (3.8) 
For c = -I we have faster convergence in Lemma f3 . 1 1 since a? = 0. 

3.2 Edgeworth Expansion of 

From the similarities between the two expressions in (|2.4p and (|2.9p . the derivation 
of -F^+i 2 follows exactly the same steps as the previous one, they would differ only 
by some constant terms. 

We see that the corresponding inner products in terms of the a's are 

(ai,/3i) = — —uqU-^ , (a5,/35) = ^ 2tf;i , (a6,/36) = 2 2m;i, 

23 23-10 23-10 

_ 2 _ 2 

77, 3 77/3 

(a2,/?2) = -2 , (a3,/53) = ^^ 2(m2 - MiWo + Mo^^i - Mo^^o + "0^0), 

23 - 10 23 - 10 

_ 2 _ 2 

77. 3 n A 2 

(a4,/54) = -2 2u2 , (07, /^t) = 1 (18 - 20c + 5c )i;o = (a8,/58)- 

23 - 10 23 - 10 



If we set 



L , ^ c-2 . 2I 
a, nis) = — -, — , O9 = — and 
c,2V ; 2^ 10 



E^2{^) = — 3^2 + (5c — 20c + 18)t;o + UiVq — UqVi + UqVq — UqWq, (3.9) 
our formula reads 

Lemma 3.2. For tc = An + 2a + 2c + 2{2n)h, 

F^+i,2{tc) = F2{s){l + a^,,uo{s)n--A+b^E^^,{s)n-l} + 0{n''). (3.10) 
For c = 2 we obtain a faster convergence as 2 = 0. 

3.3 Fine tuning 

To complete this analysis, we need to find values for the constant c = Ci in Lemma 
13. II and c = C2 in Lemma f3. 21 for which E^^ 2i^) = F^^^ 2i^)- Which is equivalent to 



-20c? + 20ci - 2 = 5c^ - 20c2 + 18 ^ (ci - -)2 + -(c2 - 2)^ = - (3.11) 
This suggest the following change of variables. 

cg = ci and cl = — . Therefore 
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dnZH) changes to 4 + 4 = - , a^^ ^ = cg , a^^,2 = '^''Cl , 
-20cl + 20ci - 2 = -204 + 3 , 5cl- 20c2 + 18 = 20^ - 2 , 



= (2(n + l) + 2cG)5 + 2-5n-ss and = 4(n + l) + 2a + 4cL + 2(2n)3s . (3.12) 
We can now give a scaling of t in terms of the size of the matrices. 

= (2(n + 1) + 2cg)^ + 2~2 (n + l)~6s for the Gaussian case, and 

= 4(r2 + 1) + 2a + Acl + 2{2{n + l))^s for the Laguerre case. 

Since all the functions derived so far are all differentiable, this new scaling will only 
change the error function but not its order and class. To keep the notations light, 
we will use the same variable to represent these error functions. We therefore have 
Theorem 11.21 and Theorem 11.31 for an ensemble of [n + 1) x [n + 1) matrices. 

Theorem 3.3. For x = {2{n + 1) + 2cg)^ + + 1)~^X and 

y = {2{n + 1) + 2cg)5 + 2-^n + 1)"^^ 

y) dx = {AAi(X, Y) - cgM{X) Ai(F)(n + 1)^^ + 

^ \x + Y) Ai'(x) Ai'(y) - + xr + f^) Ai(x) Ai(y)+ 
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zucg '^.(Ai'(x) Ai(y) + Ai(X) Ai'(F)) 



(n + i)-i + o((72 + i)-i)E(x,y)|dx 



2 

(3.13) 

The error term, E{X^Y), is again a kernel of an integral operator on L'^{J) which is 
trace class for any subset J of the reals that is bounded away from minus infinity. 

Taking the limit as F ^ X in (HrTT?!) give the one point correlation function pn+i- 

Corollary 3.4. 

For X = [2{n + cg))'^ 



22n6 



2-^n-^Pr,{x) = 2-^n-^K,,{x, x) = [Ai (X)]^ - X[Ai{X)f - cg [Ai{X)fn~^^ + 

^ {2X[Ai'(X)]2 - 3X'[Ai{X)]' + (3 - 204) Ai'(X) Ai(X)} n"! + 0(n~i)F„(X) 

(3.14) 

Note that for cg = 0, this is formula (72) of [7j. Figure illustrates the accuracy 
of equation ()3.14|) . 

For the Laguerre case, 
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Figure 3.3.1: For Cc = and n = 40, The dashed curve is the usual approximation of 
the one point correlation function from the Airy kernel (which is the first two terms 
in (|3.14j) ). the solid curve is the exact scaled one point correlation function and the 
doted curve is our approximation (|3.14|) 

Theorem 3.5. For x = A{n + l) + 2a + Acl + 2(2(n + 1))^X and 
y = 4:{n + 1) + 2a + 4:Cl + 2{2{n + 1))^Y with X,cl and Y bounded, 



The error term F{X, Y) is again the kernel of an integral operator on L'^{J) which is 
trace class for any subset J of the reals which is bounded away from minus infinity. 

Taking the limit as y goes to x in ()3.15j) gives the one point correlation function 
p'^_^_i in the Laguerre case. 



K^^,{x, y) dx = [k^,{X, Y) - 2hcL Ai(X) Ai(y)(n + 1)^^ + 



10 L 



{X^ + XF + Y^) Ai(X) Ai(F) - (X + F) Ai'(X) Ai'(F) 




Corollary 3.6. For x = 4(ra + c^) + 2a + 2(2ra)3X , 



2{2n)\pl{x) = [Ai'(X)]2 - X[Ai(X)]2 - 2^ CL[Ai{X)]'n--^ + 
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[3X2[Ai(X)]2-2X[Ai'(X)]2+2(l-10c' 



i) Ai(X) Ai{X)]n-l+0{n-')F{X) (3.16) 
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Note that for cl = —a/2 this is formula (73) in [7]. Figure 13.3.21 illustrates 
the accuracy of our result. Combining Lemma 13.11 Lemma 13.21 and the fine-tuned 
constants in (j3.12p give Theorem 11.41 




Figure 3.3.2: For = 0, a = ^ and n = 40, the solid curve is the scaled one 
point correlation function, the doted curve is our approximation from ()3.16|) . and the 
dashed one is the first order approximation from the Airy kernel. 
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A Proof of Theorem ll.il 

In this section we will follow Szego |TH], Section 8.75. With the following changes. 



• We introduce a new variable c in the definition of /„ to fine-tune the final 
result. In = (4n + 2a + 2c) 2 instead of /„ = (4n + 2a + 2)2 . 

• We will use the second order Hankel expansion of the Bessel function of large 
argument instead of the first order. 

• We define = In + (2/ri)~^t instead of i = In — {^ln)~^t to match the 
definition of the Airy function with the one commonly used. 

• We give an estimate of the error term as a function of the independent variable 
recalling the generating function of the Laguerre polynomials. 

n=0 



anrl fnr f^r\^f^ cr'a liner / = I An -U Vn/ 4- 

4 ' 

uj ^ InZ. we deduce from ()A.1|) 



The substitutions x — ^ and — ^ — and for edge scaling In = (4n + 2a + 2c) 2 



rfn + a + l) 4' 27rz 



J-y 

(A.2) 

where 7 is a symmetric contour enclosing the origin. 

Using the Hankel expansion of the Bessel function of large argument (see, for 
example, Olver [13 pgs 130-132] ): that'^ as \lnz\ — > 00 in | arg(z)| < vr — 6{< vr). 



e 2 J„(e 2 ^l^z) = i2ii^ln) 2 

Aa —1 . g;^^ ^/ 



Substituting this in ()A.2|) for ^ = /„, + (2/„) at with t bounded, 
LZie) , 1,„ 2-(/„ + (2/„)4t)-^ 



r(n + a + l)^ 4^ 27ri 



/-^"-"(27r(/„+(2/„)'^t)/nr- 



^In the following formula the O-tcrm is actually of the form 
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,,2 _ 1 ,2 .2 



A 8/2 -(- 9-Il 



^/2^(l + 2-3/„3^) 



The ± signs would be taken according to which one gives the larger contribution. 

1 _4 

To simplify notations we set p = (4a2 _ l)/(8/2^(l + 2-3/„ H)). 

The integral over 7 can be split into the path in the upper half-plane and the 
path in the lower half-plane. However the lower-half plane part of the path can be 
transformed into the upper half plane via the transformation 2; — > 2, taking into 
account the orientation change and the fact that the integrand is an analytic function 
of z on 7. The integrand will be transformed into its complex conjugate under this 
transformation. The Jacobian of the transformation is —1 so it will reverse the 
orientation once again. Thus the lower half plane contribution of the integral is equal 
to the conjugate of the upper half contribution such that the integral over 7 is exactly 
twice the real part of the integral over the upper half portion 7+ of the contour. Thus 
flA.3|) is equal to 

2"(/„ + (2/„)-^t)""C2n-a(2^(;^ ^ ^2/„)-it)/„)-| ■ 27^e{^G + + K} (A.4) 

where 

G= f (1 - p)e-''/i(-)+2"^'l-t z'-^ dz= I e-''f'^'^g{z) dz (A.5) 
"'7+ "'7+ 

z^ 1 
with f^{z) = — -z + - \og{z) 

H = e(°+5)-i /" (1 + p)e-'"/2W-2-3«|.t ^c-l f e-'"^2(^)/i(2) dz (A.6) 

2^ 1 
with f^{z) = — + z + - \og{z) 

and i^ = 0(/-^) I |e-%^-|'nM^)±'"-±2-^i|.t| |^^| ^^.7) 

We will use the steepest descent method to find an asymptotic expansion for G, H 
and K for large keeping a fixe and c bounded.® 



^Actually c can grow with n, but wc arc not going to look at this problem. 
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A.l Steepest descent method for G 

The steepest descent for this section is with respect to large In in ()A.5j) . The saddle 
point condition f (z) = gives zq = 1 with / (1) = 0. 

The steepest descent curve leaves ^tq = 1 at angles 0, n/S, or 27r/3. The direction 
of maximum decreases of /i is 27r/3. We deform 7_|_ at zq = 1 such that the resulting 
contour leaves 2:0 = 1 at angle 27r/3 as the line segment 

z = 1 + 2kn^pe^ with < p < where < 6 < - (A.8) 

6 

then along the segment symmetric to this segment with respect to the imaginary axis 
and finally connect the tip of these two line segments by an arc of circle centered at 
the origin with radius r. For simplicity, we will also call this path^ 7+. 

In the following sections, we will estimate the contribution of each portion of the 
contour to G. 

A. 1.1 On the arc of circle z = re*^* 



Since 7le{— fi{z)) = —^r^ cos 20 + r cos0 — |log(r) is a decreasing function of for 
< < TT, the major contribution of this arc is bounded above by the value of the 
integrand at the end point where p = n^. A deformation of the path of integration 
near —1 is immaterial so the estimate of the remainder of the path is again bounded 
above by the value of the integrand at p = n^. (Therefore in the next section we will 
focus on the asymptotics on the line segment in the first quadrant only.) 

2 1 _2 

If we set z = re*^ where = 1 + 23/„ ^n^^ — 2^1^ ^n^ and 9q < 9 < n — 9q 
where 9q is the angle that the ray from the origin to the tip of the line segment 

1 _2 ^ 

makes with the real axis, cos(^o) = ^"^^^'r " • modulus of this integrand in ()A.5|) 
is of order 

3 1 1 1 — — 

^c-2gfc(9) ^i^gre k{9) = --r"^ cos(26') + r cos(6') - - log(r) + 2"3/„ 3 ^ ^ cos(6') 

_1 -4 

k' has roots vr, and ^2 where cos(6'2) = ^^'^ J'" * , for large /„, < ^2 < ^0 • 
k increases from to 92 and decreases from ^2 to vr. Therefore the maximum of 
k{9) on the arc of circle is at 9o. We can also use this as an upper bound of the 
contribution of the contour around —1. Thus the contribution of the arc of circle and 
the line segment on the second quadrant is of order of the modulus of the integrand 
evaluated at the tip of the line segment in the first quadrant. We can therefore focus 
our attention in the next section on just the line segment in the first quadrant. The 
parametrization ()A.8|) shows that p will go to infinity with n, so we want to estimate 
the contribution not only of the line segment but of the whole ray < p < 00 . The 

error that we make by taking the ray is | ^'^"se"'"-^^*^^)"*"^ ^'•"'^^dz] . 

^See fHl Section 8.75] for an illustration of this path. 



20 



On the ray, -llh{z) = f/^ _ /L + 0{p^)ln~' , thus 

We therefore need to give an estimate an estimate of the integral on the right of this 
last equality and an estimate of the integrand when p = in ()A.5|) . 

Integrand when p = 

1 — - 3 2 - § 

The order of the integrand for large n is 0(1)23Z„, =^64'"+^ g-^+pie~3~_ 
We need to give an estimate of the factor u := e~"3"+''*'^^ as a function of t. 
If t is positive, it is of order e~~ ■ e 2" or of order e 3~ ■ e . 

3 

If t is negative, we have \u\ = e~~~^~ where y = —t. For trace class convergence we 
need an estimate that will decay exponentially for large t. The expansion of the expo- 
nent in \u\ around his critical point po = ~\/\ is ~\{\)^ ~^ \/%{P~^ ~\{P~^ ■ 

_{r^ + ^/^ _2 3 „3a 

Thus for p = n , u is of order e 3 ■ e~3'.2) or of order e 3 ■ since y is 



bounded and positive. So in either case the contribution is of order, e 



„3i5 



V . P-I*l 



Tail integration 

3 ^ 

We show that the contribution of j^s e~^~^dp is of the same order. 

3 3 3 

If i > , /J e-^-^dp < /J e-'r-idp = d /J e'^rdp < e'^ /J e'sdp 
so it is of order ■ e~2 

If t is negative, a similar change of variable and expansion of the integrand leads to 

/J e-4-^dp = /J e-i(i)^+v^(''+v^)'-^(''+v^)'rfp 

e-|(|)^ /J e-^(^+v4)^ni-3(p+v^)-v^]^^ ^ ,-§(|)i jj e-^(^+v^)'rfp 

This is of order, e'i^a)^ e"^^^'^v^^c/p or of order e'i^a^^ ■ e~^^'"''"'"v^^ 

So the contribution of this integral is of order e~~ ■ e"'*' 

We conclude this subsection by recording that the error that we make by neglecting 
the remainder of the contour and considering the integral from zero to infinity instead 
of zero to is at most of order 

2^Z-^et'"+2-*/lv4.e-l*l (A.9) 



2 

A. 1.2 On the ray 2 = 1 + 2h^^pe^, p G [0, 00) 
The Taylor expansion of /i at 20 = 1 is 

Mz) = -! + iln' - I ET=.c,,,ipl-h' with c^,, = i-ir-^ , and 
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2 slntz = 2 s/nt + pte^^\ Taking in account the error estimate from ()A.9|) . 
the substitution of these quantities in G give, 

2 1 2 f' /"OO 3 27ri 6 1 



gnip) = (1 ^2 -^)(l + 23/„Ve~) 2e5^fc=4Ci.*(p'" ) 

8/2(1 + 2^/„Ve^)(l + 2-^/n^t) 

(A.IO) 

A. 2 Steepest descent for H 

The analysis for H differs from that of G in the location of the saddle point, and 
the orientation of the contour. A similar analysis shows that the saddle point is now 
at zq = —1, the final contour of integration is the same but oriented in the opposite 
direction. It leaves zq at angle 7r/3. The error estimate on the arc of circle and on 
the tail of the corresponding ray is the same. The new parametrization on the ray is 

z = -l + 2k~hfp = -(1 + 2knh-^p), < p < oo 
The Taylor expansion of /2 at z = — 1 is 



f2{z) = -l + l^ + ln'J-lJ2''^,kp'ln^\ With C,,,= {-lf—^ 

k>4 

This leads to 



k 2 7vki 
3 



7 / \ — 1 1 2 27ri,. 3 1 v^oo . . ^-skj2 

K{p) = (1 -2 -^)(l + 23/„3pg-— )C-5g5l^,^4C2,fe(pi„ ) Z„ 

811(1 + 2kn'pe~'-f){l + 2-knH) 

(A.ll) 

A. 3 Asymptotics for K 

The asymptotics of the integral factor in K depends on the leading term in the 
expansion of either G or H, depending on which one is larger as shown in ()A.7|) . 
But from the previous analysis, the leading term of both G and H are of the same 
order. Thus K is also of order of a nonzero linear combination of the leading terms 
in G and in H times 0(/~^). In our case^° we take the linear combination to be 
7ee^(G + H). 
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The choice of this representation of the error is for trace class convergence of the final result. 
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A. 4 Conclusion 

Note that from ()A.10|) and (jA.llI) we see that except for the 0-term, H = —G. 
The change of variable p i— > pe~~^ transform gn{p) into a real function gi{p) and 



G = 2kn~'ei''-+^'^'^' ■ ( / e^-'^+P'^gi{p)dp + 0(e-'^)e-l*l), note also that 

Jo 

27^e— (G + H) = 2- 2h~h-^''+^~'^'"'- 
2Tii 



2Tvi 

^' e^-'T+pi^g^{p)dp+ e(-V+p*)(^,(p)rfp] + 0(e-'r)e-l*l J> (A.12) 



27ri 



At this point we can give an estimate for K based on this last formula since the 
leading term of (^i is 1. 

2 1^1/"*^ 3 /•ooe~3~ 3 

K = 0(/^) ■ 2i/;;^et'"+2""'" * ■ — ( / e(- V+^*)dp + / e^" 

2vri Jooe~^ Jo 



27rz 



q /"ooe 3 



where Ai{t) = ^— ( / "3 ~^P^^dp+ I "3 is the Airy function 



ooe 



(A.14) 

To simplify notation, we will combine the two paths of integration in G and call 
the new path a. 

With this notation, ()A.4|) gives 



e-^L^^(e') = (-l)"22"r(n + a + l)2"(0""e-^C'"~"(2<^n)~^2 ■ 2kn ^ei'^+^-nzlt. 

[1 /■ e(-4+''*)^i(p)rfp + 0(OAz(t) + 0(e-^)e-l*l] (A.15) 
27r2 

Using ()A.10|) and the help of Mathematica, we have the following expansion^^ of 

_ 2 

gi in powers of ^ • 

i.p' . 3, , ,-| 2/ p5 p8 (c-|)p5 (c-|)p2(c-|)\ _4 

1 + 23 — + c - - p „ ^ + 23 + + ^ ^ + ^ ^ „ 3 

V ^ T V 2^^^^ V5 32 4 2 / 



The expansion is valid for this derivation up to In as the error estimate for K indicates. We 

-42 

choose to stop here at In ^ note also that the last term will served to estimate the error. 



23 



, ,- (c-|)/(c-|) , p6 P'p''h 3,3, 5,, 7, 
+ ( ^ + 7-10 + 192 + 3(^-2)^ (^-2)^^-2^ 



2 8^ 2^ 5 327 " 
ol( 2p^ ^ ^ ^ ^16 ^ ^ 1 / (c-|)p3(c-|)(c-|) ^ 1\ 



7 12 25 80 3072' 4 V 3 2 

3 , 1 o , 5 , , p^ , , 3 , 1 , p^ p^ p"*^^ , 

+ c - - 23p2 c - - -— + — + c - - 23p^ - — + — 
^ 2^ ^ 2^^ 5 32^ ^ 2' ^^^3 10 192' 
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+ - \)p\c -\){c- l){c + (-^ + i)(c - ^)2ip - i(-4a^ + l)2ip ) 
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+q{p)0{ln ^ ) for some polynomial q 
The integral in ()A.15j) can be expresses as a linear combination of the Airy function 

3 

and its derivative using Ai'-'^''(t) = ^ J^G~^~^''^p'^dp and this expansion of gi. Using 
the Airy differential equation Ai"(t) = t Ai(t), it reduces to an expression involving 
only the independent variable t, Ai and Ai' . 

The contribution of the last term of the expansion of gi is a finite combination of the 
form X]m=o(^'n^(^) ^K^) Ai'(t)) for some polynomials p and g. If t is bounded 
away from minus infinity, this is of order Ai(t). In this paper we assume therefore 
that this is the case for t. 
If ^ = /„ + (2/„)~3t we have 

= ^-(-+1)^4-2-^4*^1 _ ii2i/;i + (|i - {2a + 1)^)2§/;^ + (-3^;, + 



{2a + l)i)l-' + 0{l-H% 
Stirling formula gives T{n + a + l) = 257r5n"+"+5e-"[l + ^°^"+^^)+^ + 0{n~^)] and 
we have C'""'"-^ = 2-2"-2«-f n— ""le"^ [1 + i-+c)i-9»+3c-w) ^ o(^-2)]_ 

Thus 

{-l)^2'^^r{n + a + l)2"(0""C^"""(2</„)-52 ■ 2k~h^^'^^'^^"'e-^ = 

(_l)n2"a-|^-|(l + 3a!+2a^6ac±3c^^ 

{-^ + i2a + l)i)l-' + 0il-H% 



This in ()A.15|) together with the expansion of gi give the desired result. 
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^Actually 



C = (4n + 2a + 2c)^ + + 0{n- 

23^6 



but due to the smoothness of e 2 L"(^^) the error that we make by removing the O-term is 
neghgible if we are aiming for an accuracy of order n^^. 
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For ^ = {An + 2a + 2c) 5 + and t bounded 



2- 10c + 5c2-5a 
10 ■ 2§ 



tAi(t) H ^Ai'(t)ln-3 + 

20 ■ 2i 



,,5a — Ibca + 2c'^ — 15c^ — 56c — 6 c — 1 n, . , , 

[( + -^^) Mt) 



60 



40 



(c- l)(5(c-2)c-3(2 + 5a)) , 4, , 



60 



(A.16) 



This is the desired formula for this section. Figure 1X^41 gives an illustration of our 
asymptotics. 




.06 



Figure A. 4.1: For a = — c = 1 and n = 40, the solid curve represents e 2 the 

dashed curve is the usual first order approximation of e~~L"{^'^) in term of the Airy 
function (Which is the first term approximation in ()A.16|) ). the doted curve represents 
our approximation. These are functions of t where C, = {An + 2a + 2c) 2 + 
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